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While  the  contents  of  this  report  are  considered  to  be 
correct,  they  are  subject  to  modification  upon  further  study. 
This  report  does  not  promulgate  official  Air  Force  policies 
or  positions.  The  technical  conclusions  are  solely  those  of 
the  author. 


ABSTRACT 

This  report  presents  a new  method  to  compute  a more  nearly 
optimal  initial  basic  feasible  solution  for  the  Transportation 
Model.  The  inte;?gration  of  two  techniques;  (1)  The  Decision 
Index  (DI)  and  (2)  An  Admissability  Index  (AI) , have  resulted 
in  50  to  75  percent  reductions  in  computer  run  time  required 
to  derive  an  optimal  solution. 


FOREWORD 


This  report  and  the  BEST  program  were  prepared  by  the  Modeling 
Branch  of  the  Air  Force  Military  Personnel  Center  in  response 
to  the  need  to  solve  large  transportation  and  assignment 
problems  in  the  management  of  the  approximately  one  million 
personnel  employed  by  the  Air  Force.  The  work  of  Dr.  Joe  Ward 
of  the  Air  Force  Human  Resources  Laboratory,  who  originally 
authored  The  Decision  Index,  is  acknowledged  for  its  contribution 
to  the  technique  presented.  Additionally,  Gloria  Jeaneatte 
McWilliams,  whose  thesis  at  the  University  of  Texas  indicated 
the  efficiency  of  the  Decision  Index  in  developing  approximate 
solutions  for  transportation  problems,  is  acknowledged. 


INTRODUCTION 


To  solve  a transportation  problem  you  need  an  initial  basic 
feasible  solution.  The  better  (more  nearly  optimal)  the 
starting  solution,  the  more  rapid  the  transportation  algorithm 
will  go  to  convergence.  As  so  succintly  stated  by  Hadley  (#1) 

in  1962,  " , It  is  worth  while  to  spend  some  time  finding 

a 'good'  initial  solution  because  it  can  considerably  reduce 
the  total  number  of  iterations  required  to  reach  an  optimal 
solution."  For  any  transportation  problem  there  exists  a 
continuum  of  initial  solutions  ranging  from  the  Northwest 
Corner  Rule  (probably  the  worst)  to  one  of  the  numerous  optimal 
solutions.  The  objective  of  this  report  is  to  describe  a new 
technique  which  derives  initial  solutions  which  are  very  near 
optimum,  i.e.  will  require  a reduced  n\amber  of  iterations 
to  reach  an  optimal  solution. 

McWilliams  (#2)  investigated  several  promising  methods  for 
developing  approximate  solutions  to  the  transportation  problem 
in  her  thesis  at  the  University  of  Texas.  In  brief  her 
conclusions  dismissed  the  Northwest  Corner  Rule  as  a serious 
alternative  and  found  that  the  Ward  Decision  Index  and  the 
cq,  in  both  static  and  dynamic  versions,  were  superior  to 
other  algorithms  currently  in  use,  e.g.  row/column  minima, 
cost  minima,  and  Vogel's  method.  She  further  noted  that 
Ward's  method  and  method  differ  only  by  additive  and 
multiplicative  constants  in  the  static  version  when  M=N  and 
for  this  reason  the  initial  feasible  solutions  are  identical 
for  the  two  methods.  (But  this  is  only  true  when  M=N  in  the 
static  versions).  It  should  be  noted  that  McWillicims' thesis 
used  a slight  misinterpretation  of  Ward's  DI.  The  original 
DI  computation  (#3)  involves  C.j  and  Ci. , the  sums  of  the 
column  and  i^  row  of  the  cost  tableau.  These  sums  consider 
the  quotas  in  their  computation.  McWillieims' version  of  the  DI 
computes  the  row  and  column  sums  without  quota  consideration. 

Ward  suggest  (#4)  that  this  incorrect  computation  of  the  DI 
may  not  have  significantly  degraded  the  usefulness  of  the  DI 
for  starting  a transportation  problem. 

The  original  promulgation  of  DI  was  made  by  Ward  (#3)  in 
seeJcing  a near  optimal  solution  to  the  sequential  assignment 
problem  confronting  Air  Force  assignment  counselors  in  dealing 
with  job  assignments  for  non-prior  service  personnel.  Counselors 
were  required  to  make  sequential  assignment  decisions  in  the 
absence  of  perfect  information.  As  each  airman  arrived  at 
the  assignment  station,  he  was  given  his  initial  Air  Force 
career  assignment  without  knowledge  of  the  vast  labor  pool 
yet  to  be  considered.  The  result  was  intuitively  sub-optimal 
where  each  sequential  assignment  subjectively  maximized  the 
payoff  criteria  against  whatever  happen  to  be  remaining  in  the 
job  pool.  Ward  conjectured  that  given  estimates  of  the  row 
and  column  means,  of  admittedly  imprecise  data,  that  the 
sequential  assignment  problem  could  be  solved  in  a near  optimum 
manner  employing  the  Decision  Index.  Wherein  the  Decision  Index 
is  computed  thusly: 
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DI 


ID 


MC 


N(M-l) 


ID 


- Ci.  - C.j  + C, 


) 


M = Number  of  persons  to  be  assigned. 


N = Number  of  jobs  to  be  filled, 

Cij  = Productivity  of  the  ith  person  on  the  job, 

Ci.  = Cij 

j=l 

C.j  = ■§.  Cij 
i=l 

C..  = ^ £ Cij 

i=l  j=l 


As  is  discussed  in  Ward's  paper  (#3),  this  DI  is  the  mean 
value  of  the  cost  used  to  compute  the  objective  function  for  all 
[M-l| 1 possible  assignments  of  the  i^“  man  to  the  j^^  job.  (Note: 
IM-Nj  -ipj^ig  seemingly  powerful  technique  has,  to  the  best  of  this 
authors  knowledge,  not  been  widely  implemented  either  as  an  assign- 
ment tool  or  in  conjunction  with  the  Transportation  Model  as  is 
now  proposed) . 

Early  attempts  to  implement  DIs  as  a starting  solution  for  the 
transportation/assignment  problem  were  not  successful.  These 
failures  resulted  from  the  type  of  problem  being  examined  in 
which  M,  the  number  of  personnel  to  be  assigned,  far  exceeded 
N,  the  number  of  jobs  available.  In  these  first  efforts,  it 
was  found  that  a DI  start  was  little  better,  and  in  fact  more 
costly  in  computer  time,  than  a modified  row  minima  start. 

These  early  failures  led  to  the  consideration  of  an  admissability 
index. 


METHOD 


Background:  A straight-forward  application  of  DI's  to  the 
personnel  assignment  problem  where  M> N resulted  in  solutions 
not  much  better  than  the  modified  row  minima.  This,  xn 
hindsight,  is  obvious.  If  all  contenders  were  at  least 
marginally  qualified  for  all  jobs  in  the  case  where  2000 
persons  sought  300  jobs,  then  the  first  300  contenders  would 
enter  the  solution  and  remaining  1700  would  not.  This  result 
could  be  marginally  better  than  a modified  row  minima  but 
would  consume  greater  computer  resources  in  deriving.  Thus 
the  need  for  some  method  of  selecting  the  order  of  admissability 
to  the  solution. 

The  statistic,  “TJj  , possessed  an  intuitive  appeal  and  has 
proven  a useful  admissability  index  in  application  to  personnel 
assignment  problems  investigated.  Here  represents  the  mean 
value  of  the  costs  for  the  i^^  person  for  all  real  jobs  (in 
application  a shadow  job  is  used  which  has  a quota  of  M,  the 
number  of  personnel  to  be  assigned  and  is  not  used  in 
computing  ^‘■c  and  ) . The  represents  the  standard  deviation 
of  the  costs  for  the  i^^  person  for  all  real  jobs.  The  mean 
cost  for  a person  is  an  overt  index  of  which  personnel  will 
enter  the  solution.  In  general  those  persons  with  smallest 
mean  costs  will  be  in  the  solution  and  those  with  the  greatest 
mean  cost  will  not.  Inclusion  of ^4  as  a moderator  provides 
for  earlier  introduction  to  the  solution  of  those  who  might 
be  otherwise  unattractive  but  have  a few  good  assignment 
possibilities  and  limited  alternatives. 

In  calculating  the  statistics  ^'';and  vas  well  as  the 
Decision  Index,  an  important  deviation  was  introduced.  BIG  M 
was  observed  to  possess  dramatic  influence  in  the  starting 
solution.  (BIG  M originally  generated  in  the  cost  tableau 
is  two  orders  of  magnitude  greater  than  the  largest  admissable 
real  cost.)  Early  on,  it  was  observed  that  using  the  original 
BIG  M caused  marginal  assignments  in  the  starting  solution 
which  subsequently  left  the  optimal  solution.  This  of  course 
results  from  the  use  of^iin  the  admissability  index  and 
the  follow-on  use  of  the  column  means  in  computing  the 
Decision  Index.  Reference  to  an  example  best  illustrates 
why  this  occurred. 

Given  10  person  competing  for  five  jobs  with  admiss- 
ability indices  computed  using  the  original  BIG  M, 
the  following  results: 
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Person# 

AI 

1 

200 

1000 

.20 

2 

210 

1050 

.20 

3 

1100 

4400 

.25 

IN 

4 

350 

1166 

.30 

5 

300 

750 

.40 

THIS 

.50 

7 

450 

818 

.55 

8 

460 

767 

.60 

OUT 

9 

500 

769 

.65 

10 

550 

786 

.70 

The  moderating  effect 

of 

has  effected  exactly  as 

intended  and  persons 

#3  and 

4 were  driven  upward 

on  the  admissability 

list. 

Given  further  that  there  was 

one  highly 

difficult 

job  in  the 

quota  bank 

on  which  most  personnel  had 

a BIG  M cost  and  person  #3  had  a high  cost. 

marginally 

acceptable 

. The 

resulting  large  column 

mean  for  the  difficult  job  when  used  in 

computing 

the  DI  (see  discussion  below)  caused  person  #3  to 

have  an  attractive  DI 

and  to 

enter  the 

starting  solution 

on  an  assignment  that 

later 

must  leave 

the  optimal 

solution. 

If  however 

a suitable 

substitution  were 

made  for 

BIG  M,  the 

following 

results 

Person# 

AI 

1 

200 

1000 

.20 

2 

210 

1050 

.20 

5 

300 

750 

.40 

IN 

6 

400 

800 

.50 

4 

340 

654 

.52 

7 

5TS“ 

. 55 

8 

460 

767 

.60 

9 

500 

769 

.65 

10 

550 

786 

.70 

3 

800 

1111 

.72 

Note  that  #3  leaves  the  starting  solution  and  #4 
remains  but  enters  later. 

The  above  ex^^nple,  intended  only  as  illustrative  and  probably 
not  derivable  with  real  data,  demonstrates  an  anomaly  of 
Decision  Indexing  noted  by  and  currently  under  study  by  Ward (#5) 
A rigorous  derivation  of  how  to  handle  this  problem  is  expected 
from  Ward's  study  (#5). 

In  this  work,  an  empirical  derivation  of  BIG  M * 1025  with  real 
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costs  ^ 998  and  ^ 0 was  found  to  result  in  good  starting 

solutions . 

Another  important  consideration  is  the  cost  of  the  shadow  job. 

In  early  usage  of  transportation  modeling  to  solve  personnel 
problems,  shadow  jobs  were  coated  at  a fixed  cost,  less  than 
BIG  M but  greater  than  the  greatest  admissable  actual  cost. 

The  rationale  was,  of  course,  to  assign  all  eligible 
contenders  to  actual  jobs  and  assure  that  no  BIG  M assignments 
were  in  the  optimal  solution.  When  this  shadow  cost  was  used 
in  computing  DIs,  the  shadow  DIs  did  not  discriminate  over 
contenders  since  all  DIs  in  the  shadow  job  column  were  equal 
to  zero,  i.e.jCshedow  - (Cshadow)  *]yn-  Again  the  statistic 
Pci  H — ’J 

'&ci'  (the  AI)  has  an  intuitive  appeal.  Since  there  will  be 
no  contention  concerning  people  with  small  AIs,  they  will  be 
in  the  solution  and  likewise,  the  large  AIs  will  not  be  in  the 
solution,  the  cost  index  for  shadow  jobs  must  discriminate 
about  the  threshold  just  prior  to  exhausting  the  quota  bank. 

At  this  point  there  are  perhaps  contenders  who  should  not  be 
included  in  the  initial  start  but  should  be  assigned  to  the 
shadow  job  while  the  search  continues  for  more  cost-effective 
assignments.  In  several  test  runs,  the  AI  used  as  shadow 
cost  has  demonstrated  the  required  discrimination  about  the 
threshold.  Therefore,  from  empirical  considerations  only, 
the  AI  is  proposed  for  this  role. 

Decision  Indices,  as  promulgated  by  Ward  (#3)  , were  computed 

with  Dlij  =1  (MCij  - Ci.  - C.i  + C..)  (see  page  S’ 

N (M-1) 

above  for  amplification) . Several  elements  of  this  calculation 
are  not  needed  in  computing  DI  for  starting  a Transportation 
Model.  Since  the  initial  assignments  are  made  in  a sequential 

manner,  the  division  bv  i may  be  deleted.  Also  the 

' 'T(M-l) 

subtraction  of  the  row  total,  Ci*,  may  be  removed  as  well  as 
the  addition  of  C..,  the  total  sum  of  the  cost  tableau.  This 
leaves  MCij  “ ^*3  which  is  computationally  easier  to  implement 
as  Cii  - C.i  , or-the  deviation  from  the  column  mean. 

M 

Software:  The  BEST  program  is  presented  in  ALGOL  in  Appendix  A. 

The  various  computations  executed  in  BEST  to  produce  the  basis 
for  Langley's  Primal  Simplex  Transportation  Model  (#6)  will 
be  documented  in  this  section.  As  appropriate,  reference 
will  be  made  to  line  numbers  in  the  program  listing.  In  line 
Joo  through  >{7  9 0 , declarations  of  files  and  variables 
are  found.  Only  three  of  the  files  are  pertinent  to  the 
general  application  of  starting  a transportation  problem. 

These  are  INPl  or  CAREERS /DATA/ TRANS  which  provides  the  number 
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of  destinations,  the  number  of  sources,  the  demand  at  each 
destination,  and  the  cost  tableau.  ADVAN  or  CAREERS /DATA/BAS IS 
which  is  used  to  write  the  basis  for  later  input  to  Langley's 
Primal  Simplex.  The  third  file  DUS  or  CAREERS /DATA/REVERSE 
is  used  as  a transient  storage  medium  in  computing  the  DI's 
from  the  cost  tableau. 

At  line  number  l^UOO  , the  size  of  the  problem  is  read  in  from 
INPl.  These  data  are  then  printed  in  the  output  report. 

Line  numbers  S’ooo  through  to  ioo  provide  more  declarations. 
Three  defines  are  shown  at  lojoo  to  lO  600*  These  are  used 
to  bit  pack/unpack  the  variable  A which  is  tag  sorted  in  the 
procedure  SORTER. 

An  in  place  tag  sort  is  shown  in  the  procedure  SORTER  at 
line  number  to  . 

Demands  at  each  of  the  destinations  are  read  into  the  array 
NJ(^S  at  line  number  • Then  the  total  number  of  jobs 

to^llocated  is  accumulated  in  the  variable,  TOTJOBS . 

Costs  are  read  in  sequentially  from  INPl  at  line  numbers 
2,  4 s' (TO  through  3o?o0  . jn  this  block  of  code  the  BIG  M 

adjustment  is  made  and  then  the  necessary  computations  for  com- 
puting /^ci  and  are  made.  From  these  calculations  the  AI 
is  computed,  bit  packed  with  the  row  index,  and  stored  in  the 
array  ROWTOT.  The  AI  is  also  stored  in  cost  (0)  to  be  used 
later  as  the  cost  for  the  shadow  job.  The  adjusted  costs 
are  then  written  to  the  transient  file,  DUS. 

SORTER  is  invoked  at  line  number  30^)00  to  sort  the  array  ROWTOT 
on  ascending  AI  value. 

Beginning  at  line  number  ?2.3oo  through  line  number  , the 

adjusted  costs  are  randomly  read  (based  on  the  sorted  AI)  from 
the  file  DUS,  the  DI's  are  computed,  and  the  basis  is  stored 
in  the  appropriate  arrays.  Of  particular  note  in  this  area 
is  the  decrementing  and  testing  on  the  variable  TOTJOBS  which 
is  the  number  of  jobs  to  be  allocated.  When  all  jobs  have 
been  allocated,  this  area  is  left  and  a wrapup  area  is  used 
to  assign  all  remaining  personnel  to  the  Shadow  Job.  Also 
note  at  line  number  jHooo  , the  DI  is  computed  dynamically  as 
a function  of  the  number  of  personnel  remaining  to  be  assigned. 

To  accomplish  this,  the  column  totals  are  adjusted  appropriately 
at  line  number  3 8u«o  through  "iiSao  after  each  assignment  in 
the  starting  solution. 

The  wrapup  area  is  shown  at  line  number  i4.oX.00  through  14SX0O  . 
Here  the  necessary  elements  of  the  basis  are  calculated  to  assign 
all  remaining  personnel  to  the  Shadow  Job.  Then  the  basis  is 
written  to  file,  ADVAN.  Other  reports  and  timing  statistics 
are  also  written  at  this  point. 
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RESULTS 


A typical  application  of  AI/DI  starting  technique  is  its 
production  usage  in  the  Air  Force  CAREERS  JOB  FINDER  System. 

In  this  system,  first-term  airman  who  do  not  have  reenlistment 
quotas  in  their  current  AFSCs  are  costed  against  jobs  which 
do  not  have  sufficient  applicants  to  meet  Air  Force  career 
manpower  objectives.  This  produces  a cost  tableau. 

In  the  problem  run,  1908  prospective  reenlistees  were  optimized 
against  939  jobs  in  79  job  categories  producing  an  implicit 
cost  array  of  1,791,612  cells.  The  explicit  array  contained 
150,732  cells  with  100%  density,  i.e.  inadmissable  person- 
job-matches  were  costed  at  BIG  M but  were  defined  to  Primal 
Simplex. 

Primal  Simplex  using  Langley's  internal  artifical  start 
converged  to  a solution  in  1,643  seconds  of  processor  time. 

The  problem  was  then  solved  using  the  AI/DI  to  generate  the 
initial  basic  feasible  solution.  This  required  87  seconds 
to  generate  the  basis  and  306  seconds  in  Langley's  Primal 
Simplex  to  go  to  an  optimal  solution.  That  is,  the  problem 
takes  4 times  longer  to  run  with  the  artifical  start  than 
it  does  with  the  AI/DI  start. 
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CONCLUSION 

Hadley's  truism  staiids.  Ward's  contribution  to  the  arena  of 
optimization  •' s again  recognized.  McWilliams's  initial 
investigation  has  served  to  stimulate  further  investigation. 

The  AI/DI  technique  is  implemented  with  the  prototype  code 
shown  at  Appendix  A.  Results  are  satisfactory.  Large  person- 
job-match  problems  are  now  taken  in  stride  without  regard  to 
ADP  impact.  In  fact,  these  problems  consume  so  little  time  that 
the  stringent  limitations  on  memory  utilization  can  be  relaxed. 
This  relieves  the  user  of  the  onerous  problems  of  memory 
management  which  may  prevail  in  some  installations.  There  are, 
however,  areas  in  which  further  research  may  be  fruitful. 

Given  that  many  optimization  problems  deal  with  imprecise  data, 
even  to  the  extent  of  being  subjective  in  some  cases,  how  near 
optimal  must  a solution  be,  to  be  good  enough?  Does  the  AI 
coupled  with  a dynamic  DI  get  close  enough  for  most,  some,  or 
any  optimization  problems?  Can  tests  be  applied  to  the  AI/DI 
objective  function  to  determine  if  it  is  good  enough,  i.e.  1, 
2,3,  or  perhaps  22  standard  deviations  from  the  mean  objective 
value?  Can  the  AI/DI  be  applied  to  sparse  matrix  problems? 
Network  problems?  Is  the  documented  computation  of  the  AI  the 
correct  method?  The  best?  or  only  one  of  many?  Is  there  a 
better  index  than  the  AI  to  use  for  costing  Shadow  Jobs?  Is 
there  adequate  payoff  from  the  introduction  of  Shadow  Jobs  to 
make  it  worthwhile? 
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